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Abstract 

We investigate a new type of dark matter witli couplings to ordinary matter 
naturally suppressed by at least 1 order of magnitude compared to weak interac- 
tions. Despite the extra-weak interactions massive particles of this type (XWIMPs) 
can satisfy the Wilkinson Microwave Anisotropy Probe (WMAP) relic density con- 
' straints due to coannihilation if their masses are close to that of the lightest state of 
the minimal supersymmetric standard model (MSSM). The region in the parameter 
^ . space of a suitably extended minimal supergravity (mSUGRA) model consistent 
with the WMAP3 constraints on XWIMPs is determined. Plots for sparticles' 
masses are given which can be subject to test at the Large Hadron Collider. As 
an example for an explicit model we show that such a form of dark matter can 
arise in certain Z' extensions of the MSSM. Specifically we consider an Abelian ex- 
tension with spontaneous gauge symmetry breaking via Fayet-Iliopoulos D-terms 
in the hidden sector. The LSP of the full model arises from the extra U{l)x sec- 
tor with extra-weak couplings to Standard Model particles due to experimental 
constraints. With R-parity conservation the new XWIMP is a candidate for cold 
dark matter. In a certain limit the model reduces to the Stueckelberg extension of 
the MSSM without a Higgs mechanism, and wider ranges of models with similar 
characteristics are easy to construct. 
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1 Introduction 

The nature of dark matter [Ij and dark energy continues to be one of the primary 
open questions in both particle theory and cosmology. It is now widely believed 
that dark matter must be constituted of particles outside the standard list of known 
particles. Chief among these are the so called weakly interacting massive particles 
(WIMPS). Supersymmetry with R-parity conservation leads naturally to such a 
particle in the form of the lightest supersymmetric particle (LSP). In the frame- 
work of SUGRA unified models the lightest neutralino is a particularly attractive 
possibility. 

In this paper we investigate a new possibility in supersymmetric models, where 
the dark matter candidate has extra-weak interactions with matter; interactions 
weaker than weak interactions by at least one order of magnitude. We refer to 
these particles as XWIMPs. We will show that such a possibility can occur natu- 
rally in certain extensions of the minimal supersymmetric Standard Model (MSSM) 
by Abelian gauge symmetries V(V)x- One such class has been analyzed recently 
[21 [31 m [5l [HI [7] . Models of this type are based on the Stueckelberg mechanism 
which arises quite naturally in string and D brane models. Further, some of the 
specific features of the models of [21 [3l HI [5l El [7] and of the type discussed here 
may have a string realization [8]. A more detailed discussion of the motivation for 
such models may be found in the above references. But there may be a wider range 
of models where extra-weak dark matter can appear. The relic density analysis of 
XWIMPS requires careful study which will be discussed later in this paper. 
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The basic elements of the models we discuss are exhibited in Fig.([T]) and involve 
three sectors: [i) a visible sector where the fields of the SM or the MSSM reside, 
(n) a hidden sector which is neutral under the SM gauge group, and {in) a third 
sector [9] which is non-trivial under the SM and the hidden sector gauge symmetry. 
Aside from gravity, the fields of the visible sector and the hidden sector commu- 
nicate only via this third sector which we therefore call the "connector sector". 
Interactions with hidden particles can of course modify predictions of the SM and 
are thus highly constrained by the precision data from LEP and Tevatron, see e.g. 

In the following, we construct explicit simple models where the gauge group 
in the hidden sector is just an Abelian U{l)x with spontaneous breaking and 
a massive Z' gauge boson. Such Abelian extensions of the MSSM occur in a 
wide class of models including grand unified models, string and brane models 
[Tot [TTl [T2| [HI [131 [HI [IS] • The explicit elements of our first example are as follows: 

1. The visible sector contains gauge, matter and Higgs superfields of the MSSM 
charged under the gauge group SU{3)c x SU (2) l x U{1)y, but neutral under 
f/(l)x. 

2. The hidden sector contains the gauge superfield for U{l)x which are neutral 
under the Standard Model gauge group. 

3. The connector sector contains the chiral fields 0"*^ with charges ±Qx under 
U{l)x and iY^ under f/(l)y. They thus carry dual quantum numbers. The 
fields in the visible and in the hidden sectors can communicate only via 
couplings with these connector fields. 

Spontaneous breaking of the U generates a mixing between the hidden and 
the visible fields. We will implement this breaking via Fayet-Illiopoulos D-terms 
[TB] . The parameters that measure the mixing are highly suppressed because of the 
precision constraints on the electroweak predictions. Their smallness is responsible 
for the extra-weak interactions of the hidden and the connector fields with the 
fields in the MSSM. 

If the LSP of the U{l)x sector which we call the XLSP is lighter than the 
LSP of the visible sector, it will be the LSP of the whole systemQ With R-parity 
conservation, it is then an XWIMP candidate for cold dark matter. 

We will also show that the above class of models reduces in a certain limit 
to the Stueckelberg extension of MSSM introduced in |3j. It is also interesting to 
investigate if an XWIMP can arise in models where mixing between the visible and 
the hidden sector occurs in the gauge kinetic energy [TH] . In a supersymmetrization 
of such a model with off-diagonal kinetic terms one as well finds a mixing between 

*The alternative acronym possibility EWIMP has already been used for a different type of 
WIMP[ir. 
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the neutral fermions: i.e., the gauginos and the chiral fermions of the visible and 
the hidden sector. Thus, these models provide another class with potential of an 
XWIMP responsible for dark matter. More generally, there may be a much wider 
range of models with similar properties. 

The outline of the rest of the paper is as follows: In Sec. 2 we work out a 
U{l)x extension of the MSSM with symmetry breaking via Fayet-Illiopoulos D- 
terms mixing between the U{l)x and the U{1)y fields. The neutralino sector 
of this system has a 6 x 6 mass matrix. It can lead to an LSP with extra-weak 
interactions composed mostly of neutral fermions in the hidden sector. The model 
reduces in a certain limit to the Stueckelberg extension of MSSM (the StMSSM) 
[5]. We briefly discuss the electroweak constraints on the parameters of the model. 
As an alternative we next consider a mixing between the visible and hidden sectors 
originating from the gauge kinetic energy which works very similarly. In Sec. 3 
we analyze the relic density of XWIMPs and show that it is possible to satisfy the 
experimental constraints via the process of coannihilation. A detailed numerical 
analysis shows that XWIMPs are candidates for cold dark matter consistent with 
the recent WMAP data. The sensitivity of the analysis on the errors in the top 
quark mass under the constraints of radiative breaking of the electroweak symmetry 
is also discussed. Conclusions are given in Sec. 4. 

2 Extra-weak dark matter in Z' models 

To start with, we introduce a class of extensions of the MSSM where a natural 
mixing of the neutral MSSM fields with fields from the hidden sector appears via 
off-diagonal mass matrices. Towards the end of this section we also discuss other 
possibilities to facilitate a mixing of a very similar type. 

2.1 Broken U{l)x with Fayet-Iliopoulos terms 

A U{l)x extension of the MSSM with a Fayet-Iliopoulos (FI) D-term can lead 
in a natural manner to extra-weakly interacting dark matter constrained by LEP 
and Tevatron data. Specific U{1) extensions have been studied quite extensively 
in the literature [TUl IT^ . The features of our model were already explained in 
the introduction. The full gauge symmetry of the model is SU{3)c x SU{2)l x 
U{1)y X U{l)x- It differs from previous formulations in that a FI D-term breaks 
the extra U{1) gauge symmetry instead of an F-term. The Abelian vector fields 
consist of the f/(l)y vector multiplet (5^, A^, Db) and the U{l)x vector multiplet 
(C^, Ac, Dc) with gauge kinetic Lagrangian given bj|^ 

>Cgkin = -\b,uB^'' - i\B(y^d^\B + \dI - ^c,,c^'' - tXccr^d.Xc + \dI . (1) 

^We only use global supersymmetry here, not supergravity, and write everything in Wess- 
Zumino gauge. 
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The superfields with components ((^^, F^) are described by 



+9xDci 



if^a'^Dj^ - V2 {igxQx<P^r\c + igYY^<P+f+XB + h.c. ) 



gyD 



where F is set to zero and the covariant derivatives of the scalars are 

D^<P^ = {d^ ± igxQxC^ ± igYY^B^)<p^ . 
Next we add to the mix the FI terms 

-^Fi = ixDc + ^yDb ■ 
Ehmination of the D terms gives us the scalar potential^ 

2 



V = ^f{Qx\<l>^\'-Qx\<P 

Minimization of the potential leads to 



Y 



(2) 



(3) 



(4) 



(5) 



-) = o, 



Ig l^xQx + gUvY^ 
glQl + g^Y^ 



(6) 



We consider the bosonic sector first. Spontaneous breaking of the electroweak 
symmetry gives rise to the mixing of the neutral gauge fields (C^, S^, A'^^), with 
A^^"- (a = 1,2,3) for the SU{2)l gauge fields. In this basis the mass matrix in the 
neutral sector is of the form 









M1M2 Ml + Wgl -Wg2gY 



M1M2 

■\v'^g2gY 



ygi 



(7) 



The parameters Mi, M2 are defined so that 

Ml = V2gxQx{<p-) , M2 = V2gYY^{<p-) . (8) 

There is a single massless mode, the photon, and two massive modes the Z and 
Z' . Since ((/>"'") = 0, the superfield $^ does not enter in the mixings in the mass 
matrix for the fields in the hidden sector and the fields in the visible sector, and 
we do not consider it further. The CP-even component of the complex scalar 0~ 
mixes with the two CP-even Higgs fields of MSSM producing a 3 x 3 mass matrix 
similar to the analysis given in Ref . [3l H] . 



^As was pointed out in ^20^ this kind of spectrum together with FI couplings leads to anomalies 
in supergravity which necessitates field-dependent FI terms. We will here ignore the issue and 
only deal with global supersymmetry explicitly. 
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In the neutral fermionic sector there are two additional mass eigenstates beyond 
the four neutral fermionic states in the MSSM, Ay, A3, hi, /i2- One can reorganize 
the Weyl spinors in terms of four- component Majorana spinors xs (out of /~) and 
Ax (out of Ac) in a standard way. The 6x6 neutralino mass matrix in the basis 
iiXs, Ax); (Ay, A3, hi, 112)) reads 






Ml 


M2 











Ml 


rhx 














M2 





nil 





-Cf3SwMQ 


spSwMo 











ni2 


c^cwMo 


-s^cwMo 










cpCwMo 





-fi 








spSwMo 


-SfjCwMo 


-/^ 






(9) 



A few explanations are in order: rhx arises from the soft mass term — l^rixAxAx, 
Mo is the Z boson mass at the tree level, cw = cos9w, sw = sinOw with 9w the 
weak angle, similarly C/3 = cos/3, S/3 = sin/5, with tan/3 = {H2)/{Hi), and finally 
/i the Higgs mixing parameter of the MSSM. The mass eigenstates of the system 
are defined as the following six Majorana states 

where Xa = 1; 2, 3, 4) are essentially the four neutralino states of the MSSM and 
the (a = 1,2) the two additional states composed mostly of the new neutral 
fermions. 



We will discuss in a moment that the current electroweak data puts a stringent 
bound on e = M2/M1 such that |e| <^ 1 [6]. In this limit the masses of ^1,^2 



mgo ~ y + -rfij^ - -rhx, m^o ~ y + ^'^x + if^x ■ (U) 

For the case when the lightest of the MSSM neutralinos x° = Xi is also lighter 
than ^ = nothing much changes compared to the pure MSSM. The LSP of 
the MSSM will still be the LSP of the full system, and the dark matter candidate 
will be essentially the same as in the MSSM with minor modifications. However, a 
very different scenario emerges if is lighter than x^ aiid becomes the LSP. The 
upper bound on e translates to a suppression of the couplings of to MSSM fields 
relative to the couplings of x^ by a factor of e. Roughly speaking one can treat 
as a standard LSP x° but with couplings appropriately suppressed by at least an 
order of magnitude. This is why we call C,^ extra-weakly interacting, an XWIMP. 



2.2 Stueckelberg reduction of U{l)x extension 

In a certain limit the model of the previous subsection reduces to the Stueckel- 
berg extension of MSSM proposed in [3l H]. To achieve the reduction we assume as 
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is conventional in the analysis of MSSM that is negligible. We consider now the 
limit —>■ oo, gxQx 0, and gyY^ 0, with Mi and M2 fixed. This leads to 

^\D^<p-\' = ^{M^C, + M^B^ + d.af + ^{d.pf . (12) 

where (j)' = p + ia. The Lagrangian can be written = Cst + '^<s>+ , where is 
now completely decoupled from the vector multiplet and Cst can be written 

Cst = -\{M,C, + M2B^ + d^af-^{d,pf-zxa^d,x (13) 
+p{MiDc + M2DB) + [x{Mi\c + MsAb) + h.c] . 

This arises from the following density in superfield notation 

Cst = j (fe(fe {MiC + M2B + s + s)\ (i4) 

where C and B are gauge supermultiplets and S a chiral supermultiplet. See [3111] 
for more details. With the above one then has exactly the Stueckelberg extension 
of the MSSM. 



2.3 Electroweak constraints on mixing parameters 

The U{l)x extension of MSSM of the type discussed in Sees. (2.1) and (2.2) has 
two mass parameters, Mi and M2, which can affect electroweak physics. However, 
since the Standard Model is already in excellent agreement with LEP and Teva- 
tron data the above extensions are severely constrained. New physics can only 
be accommodated within the corridor of error bars consistent with precision mea- 
surements. For example, the SM prediction relating the W and Z masses (in the 
on-shell scheme) is given by [2T] 



Ml^ ^GpM^il-Ar) ^' ^^^^ 

In the above a is the fine structure constant at the scale = 0, Gp the Fermi 
constant, Ar the radiative correction such that Ar = 0.0363 ± 0.0019 [22], where 
the error in Ar comes from the error in the top quark mass and in a(M|). Using the 
current value of the W mass Mw = (80.425 ±0.034) GeV [22] one finds the central 
value of Mz in excellent agreement with the current data, Mz = (91.1876±0.0021) 
GeV. But the error of theoretical prediction is SMz ~ 30 MeV. Using techniques 
similar to the ones used in constraining extra dimensions [23] one may equate this 
error corridor in Mz to the shift in Mz due to its coupling to the Stueckelberg 
sector of the extended model. This constrains e = M2/M1 to lie in the range [6] 

^0.061^/1-^. (16) 



7 



A more detailed analysis of all the relevant precision electroweak parameters can 
be found in P, [7j. For our current purposes it is sufficient to know that a mixing 
parameter e in the range 0.1 — 0.01 or smaller is in principle imaginable. This 
sets the suppression factor for the couplings of XWIMPs relative to WIMPs in our 
models. 

2.4 Abelian extension with off-diagonal kinetic terms 

There is a well known example of an Abelian extension of the SM with a mixing 
between the visible and the hidden sector arising from an off-diagonal kinetic energy 
[IB]. The hidden sector in this model is called the shadow sector, the extra gauge 
factor denoted U{l)s- Specifically we write for the action £ = £sm + ^-C, where 

AC = —C^'-'C^, - ^-B'^'C^, - \D,^\' - V{(P, 0sm) . (17) 

Here is gauge field for the U{l)s, 4> is the Higgs charged under U{l)s giving 
mass to C^, and 0sm is the Standard Model Higgs. The kinetic energy of Eq. ffTTl) 
can be diagonalized by the transformation 

where a = 1/(1 - 6^y/\ ss = 5/(1 - S^)^^^- As in the analysis of Refs. H El [7] 
the mixing parameter 6 is small p4l [25] . After spontaneous breaking this type of 
model also leads to a massless photon, and two massive vector boson modes. 

To supersymmetrize the model we write the Lagrangian for the extended theory 
= -Cmssm + AC In the pure gauge sector of the theory one has 

-^-C^^'B^, - t6{XcCT^d,XB + XBCT^'d^Xc) + 6DbDc . (19) 

One can give a mass to the C^j by a Stueckelberg mechanism without mixing with 
the hypercharge as in the analysis of Ref.[5]. Thus we add a term 

A£st = j de^d9^{MC + S + Sf , (20) 

where C is the gauge multiplet for the extra U{l)s and S a chiral superfield. Ev- 
erything works very much the same way as in the standard Stueckelberg extension. 
After spontaneous breaking of the electroweak symmetry the neutralino mass ma- 
trix in the basis {{ips, X'x)', (Ay, A3, hi, h2)), obtained after rotating the Majorana 
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fermions by the use of (|T8l) . is 



Mcs 








ssCpSwMQ -Ssspsv/Ma 



.(21) 



ssCpSwMQ -c^swMq cpCwMo 
-ssS/sSwMq spSwMo -SfsCwMo 



-fx 



-fi 



The structure of the neutrahno mass matrix in Eq. (pT]) is significantly different 
from that of Eq.Q. Similar to the analysis of Sec. (2.1), in the limit ss ^ the 
states ips and decouple from the rest of the neutralinos. As before we label 
these two ^1,^2 with masses given by 



Diagonalizing Eq. fl2T]) one obtains six mass eigenstates ((■Ci; ^2)! (Xi? X2? X45 X4))- 
The situation is very similar to the models discussed in previous subsections with 
off-diagonal mass matrix. Thus we can summarize that the supersymmetrized 
model with kinetic energy mixing can also lead to an XWIMP that becomes the 
XLSP with extra-weak coupling to the Standard Model. 

From now on we use a unified notation labeling the extra-weakly interacting 
particle as an arbitrary XWIMP denoting any class of model. The small mixing 
parameter will be called e in any case and the analysis of relic density given below 
applies to all such models with XWIMPs. 



Since the interactions of XWIMPs with matter are extra-weak the annihilation 
of XWIMPs in general is much less efficient in the early universe. Thus it re- 
quires some care to ascertain if a reduction of the primordial density is possible 
in sufficient amounts to satisfy the current relic density constraints. However, the 
condition of thermal equilibrium are still satisfied for XWIMPs as long as their 
interactions are only suppressed by few orders of magnitude, say one or two. This 
requires that interaction rate P is greater than the expansion rate of the universe, 
r > H with H = T^/Mpi. For the system at hand, consisting of weakly and 
extra-weakly interacting massive particles (WIMPS and XWIMPs) the condition 
of thermal equilibrium is indeed satisfied. The XWIMPs will only slightly earlier 
fall out of equilibrium but both types of species will be produced thermally after 
the Big Bang or after inflation. This is in contrast to models where the couplings of 
dark matter candidates are only of gravitational strength or suppressed in similar 
ways. 




(22) 



3 Dark matter from XWIMPs 
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While the annihilation of XWIMPs alone cannot be sufficient to deplete their 
density efficiently such reductions may be possible with coannihilation [26]. In 
general, coannihilation could involve all the neutralinos as well as squarks and 
sleptons in processes of the type 

X'+X'j ff,WW,ZZ,Wh,--- , (23) 

where x'i = iCayXa)- Let us explain how this can potentially lead to sufficient 
annihilation of XWIMPs. 



The analysis of relic density involves the total number density of neutralinos 
n = which is governed by the Boltzman equation 



dn 



— = -3Hn - ^{cyijv){n.inj - nfn^) 



(24) 



where cxij is the cross-section for annihilation of particle species and the 
number density of Xi in thermal equilibrium. The approximation n^/n = 7f^/7f^ 
gives the well known 



where 



at 



(25) 



(26) 



the 7j are the Boltzman suppression factors 



n. 



eq 



9i 



'1 + Ai)3/2e-^- 



(27) 



Here gi are the degrees of freedom of Xii ^ = fni/T with T the photon temperature 
and Aj = (m^ — mi) /mi, mi defined as the mass of the XWIMP which is the LSP. 
The freeze-out temperature is given by 



X 



f 



In 



45 



G 



N 



(28) 



Now Nf is the number of degrees of freedom at freeze-out and is Newton's 
constant. The relic abundance of XWIMPs at current temperatures is finally 



1.07 X lO^GeV 



-1 



(29) 



Here Xf = mi/Tf, Tf is the freeze-out temperature, Mpi = 1.2 x 10^^ GeV and h 
the present day value of the Hubble parameter in the units of 100 km • s~^ • Mpc~^. 
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3.1 Relic density analysis for XWIMPs 

After all these preliminaries let us come to the specific treatment of XWIMPs. The 
naive expectation is that XWIMPs would not be able to annihilate in sufficient 
numbers to satisfy the current relic density constraints. An exception to this 
expectation is the situation of coannihilation ^26j that can drastically change the 
picture. It can contribute in a very significant way to the annihilation process. Let 
us consider the coannihilation of a XWIMP and a WIMP via the following 
set of processes 

e+e - X, 
e+x' - X', 

/ + / _ X" , (30) 

where {X} etc denote the Standard Model final states. The effective cross section 
in this case is 

^eff = Im^ ^Q + ?^)' ' (31) 

where 

Q = ^(l + A)ie-^'/^ . (32) 

Here g is the degeneracy for the corresponding particle and A = {m^o —m^o)/m^o. 
For the case at hand, the ratio cr^o^o/cr^Oj^o ~ C^(e^) ^ 1- Thus if the mass gap 
between and is large so that XfA ^ 1, then cXefr is much smaller than the 
typical WIMP cross-section and one cannot annihilate the XWIMPs in an efficient 
manner to satisfy the relic density constraints. 

If the mass gap between the XWIMP and WIMP is small and the XWIMP is 
still lighter than the WIMP we have the case of coannihilation. Let us look at a 
parameter choice with Q ~ 0(1) and Q ^ cr^o^o/cr^Oj^o. We can write cXeff in the 
form 

(Teff - a^o^o (y^) • (33) 

The result of Eq. fl33|) is easily extended under the same approximations including 
coannihilations involving additional MSSM channels. Now Eg. (1331) is modified so 



that (J^o^o is replaced by aes{MSSM) and Q is defined so that Q = Qi where 
Qi = {9i/9i){^ + Aif/^e-^f^\ When Q ~ 0(1) the XWIMP relic density is just a 
modification of the WIMP relic density modified only by the multiplicative factor 
(t^)^- It is then possible to satisfy the relic density constraints much in the same 
way as one does for the LSP of MSSM. 
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Nevertheless, the couphngs of ^ with quarks and leptons are suppressed by 
a factor of e. Thus cross-sections for the direct detection of dark matter will be 
suppressed by powers of the mixing parameter, making the direct detection of the 
extra- weak dark matter more difficult. However, will do as well as x° for the 
seeding of the galaxies. 

3.2 WMAP constraints on XWIMPs 

Extensive analyses of the relic density for WIMPS in mSUGRA|27j and its exten- 
sions exist in the literature. Recent works f2E[ [2S1 ED] have shown that the WIMP 
density in the MSSM and in SUGRA models can lie within the range consistent 
with the WMAP data. The three year data gives for the relic density [3T] 

^CBuh^ = 0.1045lHo95 (WMAP3) . (34) 
We will impose this constraint using a la corridor. 

The specific framework we consider is a Abelian extension of mSUGRA with 
a U{l)x- This means, in the MSSM we use the mSUGRA framework with the 
minimal set of characteristic parameters for the soft breaking, i.e. the universal 
scalar mass mo, the universal gaugino mass mi/2, the universal trilinear coupling 
Aq, tan/3 and sign(/i). This is our extended mSUGRA model. Its parameter 
space is subject to several constraints including electroweak symmetry breaking by 
renormalization group effects, and LEP and Tevatron conditions on the sparticle 
spectrum. The most severe of the collider data are from LEP referring to the 
light chargino mass that is expected to be greater than 103.5 GeV [32]. We 
also exhibit the bound on the Higgs mass which would eliminate rrih < 114.4 
GeV pSl although this is strictly valid only for the Standard Model. Another 
stringent constraint on the parameter space arises from the experimental limits on 
the process b s'j. The current experimental average value for the BR{b 57) 
as given by the Heavy Flavor Averaging Group [35l [36] is 

BR(6 ^ 57) = (355 ± 24t?o ± 3) x 10"^ . (35) 

In our numerical analysis we take a fairly wide error corridor and apply the con- 
straints 

2.65 X 10"^ < BR{b 57) < 4.45 x 10"^ . (36) 

For the calculation of the rehc density of XWIMPS we use micrOMEGAS 2.0[37j, 
and for the RGE computations Suspect 2.3 package of Ref. [38]. As a check on 
the results we employ DarkSUSY 4.1 [39] along with the Isajet 7.69 package of 
Ref. [10] and ISASUGRA 7 . 69. In our scans of the parameter space we implement a 
Monte Carlo technique for each scan, sampling up to 10^ points per scan. As input 
GUT scale parameters we take mi/2 G (0,1.5) TeV, mo G (0,3.5) TeV, Aq = 0, 
and multiple values of tan/? G (10 — 50), while sign(/i) is always positive. Such 
regions of the mSUGRA parameter space are within the reach of the LHC ^Ti |12] 
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and the Tevatron [13] . The particular values of tan (3 favor such a scenario. In 
the analysis of the relic density of XWIMPs we have chosen A to lie in the range 
(0,0.1). Further, the contribution of the Z' pole is included using the method of 
Ref. [Hj. However, its contribution turns out to be rather small due to the small 
width of the Z' and thus no appreciable effects occur. For the top mass the central 
value of the most recent evaluation from the CDF and DO collaborations is [45j . 

mi = 171.4 ±2.1 GeV . (37) 

The bottom mass is taken to be mh{mh) = 4.23 GeV. The allowed parameter space 
of mSUGRA is very sensitive to the assumed value of the top quark mass and of 
the bottom quark mass, and thus the dark matter analyses are also very sensitive 
to these inputs [16]. We will discuss this issue in further detail at the end of this 
section. 

In the calculation of the relic density, we find in general good agreement be- 
tween DarkSUSY 4.1 and micrOMEGAS 2.0 (up to about 15%) for values of tan/3 
in the range (10 — 40). The main result is that the WMAP3 constraints are satisfied 
by XWIMPs for a wide range of tan (3, even though the allowed parameter space 
consistent with all the constraints does depend on the value of tan (3. We consider 
two representative values in this paper, namely tan (3 = 30, 50. 

We now discuss the details of the analysis. In Fig. ([2]) we display the relic density 
constraints on the XWIMPs in the tjiq — mi/2 plane for the cases tan/5 = 30,50 
consistent with Eq. fl34p . The black region satisfies the relic density constraints 
which lie within la corridor of the central value of Eq. (1341) . while the shaded re- 
gions are eliminated due to other constraints. The other constraints arise mainly 
from the lower limit on the chargino mass and the 6 — > 57 branching ratio. The 
bound on the Higgs mass is also shown but only a small additional region of the 
parameter space is eliminated. The analysis of Fig. ([2]) shows that the relic den- 
sity is satisfied in both a low mo region, where one has typically coannihilation 
between the lightest neutralino of the MSSM and the stau, and a high mo region, 
which is characteristically the hyperbolic branch (HB) of radiative breaking of the 
electroweak symmetry , where the LSP and the next to lowest supersymmetric 
particle (NLSP) become degenerate and are mostly higgsino like. 

In Fig. ([3]) we show the parameter space in the m^o — mi/2 and the m^+ — m^o 
plane. These plots display the regions where scaling holds or breaks down which 
are also good indicators of the gaugino vs. higgsino composition of x° (the LSP 
of the MSSM). Thus in the m^o — mi/2 plot points on the straight line bound- 
ary satisfy the scaling phenomenon, where m^o ~ 0.5mi/2- Here m^o is mostly a 
Bino. More generally, scaling [18] gives rrig : m^+ : m^o ~ (6 — 7) : 2 : 1. On 
the other hand, when m^o/mi/2 is significantly smaller than 0.5 x° has a large 
higgsino component and typically arises from the HB. A similar situation arises in 
the m^+ — m^o plane. The points on the upper straight line satisfy scaling, while 
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those on the lower curved area have a large higgsino component and thus violate 
scaling. 

In Fig.(jl]) we exhibit the allowed parameter space in the rrig — plane. On 
the lower straight line along the diagonal x° is Bino-like and the scaling relation 
rrig : m^o = (6 — 7) : 1 is satisfied. Above this region has a significant higgsino 
component, and scaling is violated. Further, in Fig.(jl]) one can find the allowed 
region in the m^^ — rrig plane and, finally, in Fig. ([5]) the allowed region in the 
m^^ — and in the m/j — plot. All figures show that the permissable mass 
range for the light stop ti is rather wide, while for the Higgs there is a narrow 
window. Typically, its mass has to lie within the corridor from the lower limit of 
114 GeV up to about 125 GeV. 

Let us add a comment regarding the impact of experimental error bars on 
the top mass under the constraints of the electroweak symmetry breaking. As 
indicated above, the region in the parameter space of mSUGRA consistent with 
electroweak symmetry breaking depends very sensitively on the mass of the top 
quark, a phenomena which has been known for some time and which affects the 
relic density [4^ [5U] . We emphasize that the sensitivity of the relic density arises 
because the sparticle spectrum in SUGRA unified models, where the sparticle 
spectrum arises as a consequence of radiative breaking of the electroweak symmetry 
(REWSB) , is very sensitive to the top mass. This can be seen, for example, in the 
first paper of Ref.[4Sl| where it is shown that the stop mass can turn tachyonic 
with variations in the top mass under constraints of REWSB. However, in MSSM 
scenarios where one can fix the sparticle spectrum and vary the top mass, the relic 
density is not sensitive to variations in the top mass. In contrast, in the current 
analysis the sensitivities to the top mass arise since we are using the framework 
of SUGRA unification where the spectrum is computed via REWSB. The recent 
more accurate determinations of the top mass have now very much reduced the 
error. In the analysis of Figs. (jJHSj) we have used the central value of Eq. (l37|) . We 
now consider a la variation around this central value. Thus the results displayed 
in Figs. ([6]) are stated for rrit = 169.3 GeV, a la downward shift on the central 
value, while those of Figs. ([7]) are for rrit = 173.5 GeV which involve a la upward 
shift. Quite remarkably one finds that even a la variation with reduced error bars 
generates very significant changes in the relic density. Specifically, a lower top mass 
implies a larger portion in parameter space consistent with the constraints. 

4 Conclusion 

We have introduced a new dark matter candidate whose interactions with Stan- 
dard Model matter are extra-weak, weaker than weak interactions by at least one 
order of magnitude. Extra-weakly interacting particles can arise in a wide range of 
models, Z' extensions of the MSSM with Higgs sectors, in the Stueckelberg exten- 
sion, in extensions of the MSSM with off-diagonal gauge boson kinetic terms, and 
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possibly many other realization of small mixing between visible and hidden sector 
fields. The new XWIMPs are good candidates for dark matter if they become the 
LSP of the full system, in spite of the extra-weak interactions with the MSSM. 
They can satisfy the rehc density constraints consistent with the WMAP data via 
coannihilation. A direct observation of XWIMPs in dark matter detectors will be 
more difficult. However, indirect tests of the model are possible and should be 
investigated. 
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5 Figures 
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Figure 1: The generation of the extra- weakly interacting massive particles. An 
XWIMP is a linear combination of fields in the hidden sector and in the connector 
sector, and its interactions with the MSSM particles are suppressed. 



''Higher resolution eps figures of larger file sizes are available. 
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Figure 2: The allowed parameter space in the tuq — mi/2 plane, under the la 
WMAP3 constraint of Eq.([3l]) in extended mSUGRA for the case Aq = 0, tan P = 
(30,50) (upper, lower), sign(/i > 0), rrit = 171.4 GeV, mi/2 G (0, 1.5)TeV and 
rriQ G (0,3.5)TeV, and A in the range (0.0, 0.1). Regions eliminated by the light 
chargino mass constraint, by the light Higgs mass constraint, and by the b —>■ s + j 
constraint are also exhibited. 
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Figure 3: The allowed parameter space in the m^o —mi/2 plane and in the —Tn^o 
plane, under the la WMAP3 constraint of Eq. flS^ in extended mSUGRA for the 
same data set as in Fig. ([2]). These plots exhibit the phenomenon of scaling and its 
breakdown and in particular for the ratio m^+/m^o as discussed in the text. 
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Figure 4: The allowed parameter space in the rrig — plane and in the m^^ — rrig 
plane, under the la WMAP3 constraint of Eg. (1341) in extended mSUGRA for the 
same data set as in Fig. ([2]). 
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Figure 5: The allowed parameter space in the m^^ ~^x+ pl^iiis and in the m/j— 
plane, under the la WMAP3 constraint of Eq. llMl) in extended mSUGRA for the 
same data set as in Fig. ([2]). 
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Figure 6: An analysis for the case when = 169.3 GeV which is la below the 
central value, and tan /3 = 50 when all other parameters and constraints are the 
same as in Fig.([2l). 
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Figure 7: An analysis for the case when rrit = 173.5 GeV which is la above the 
central value, and tan /3 = 50 when all other parameters and constraints are the 
same as in Fig. ([2]). 
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